Precision Test of Quark Mass Textures: A 
Model Independent Approach 

F. Caravaglios a , P. Roudeau 6 , A. Stocchi 6 ' c 

a Dipartimento di Fisica, Universita di Milano, and INFN sezione di Milano 
b Laboratoire de 1' Accelerateur Lineaire IN2P3-CNRS et Universite de Paris Sud 
c CERN, CH-1211 Geneva 23, Switzerland 

February 6, 2002 

Abstract 

Using a Monte Carlo method, we have directly extracted from the 
available measurements, the hierarchies among the different elements 
of the quark mass matrices. To do that, we have first introduced a 
model independent parameterization for two generic class of models: 
those based on Abelian symmetries and those inspired by a U(2) hori- 
zontal symmetry. So, matrix entries are proportional to some e*, with 
e << 1 and the i's are different free exponents that we determine from 
the data through a statistically well defined procedure. We have found 
that the experimental data poorly constrain the Abelian scenarios. 

Instead, in non Abelian scenarios, these i-exponents are strongly 
constrained by the present data. We have found that contrary to a 
naive U(2) horizontal symmetry expectation, quark mass matrices 
turn out to be not symmetric. Two solutions emerge: one with 
M down > M 2 d 3 own and M 2 u f > ; and a second one with slight 

asymmetries only in the light quark sector, namely M^f < M" 2 P and 

M down > M down _ 
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1 Introduction 



Standard Model predictions have been tested so far with very high accuracy 
(at the permill level), and the agreement with the experimental measure- 
ments is very remarkable. However the Standard Model is unsatisfactory 
since all fermion masses and mixings are free and unpredictable parameters. 
Quark masses vary in a range of more than four order of magnitudes and the 
lack for an explanation is deceiving. 

If the Higgs boson is a SU{2) doublet, the four real components 

can be transformed by an SO (4) = SU(2)\ e{t x SU(2) right symmetry. The 
unique SU(2) invariant $ f $ = 0^ + 02 + 03 + 04 is also an 50(4) in- 
variant. The Higgs potential must be a SO (4) invariant, and when the 
SU (2) x U(l) symmetry is broken into U(l) em , the SO (4) symmetry is broken 
into SU(2)v = 5[/(2) le f t+right . The Goldstone bosons of the Higgs doublet 
transform as a triplet under the unbroken SU(2)y and this directly implies a 
fundamental relation between the mixing angle of the neutral gauge bosons 
and the charged/neutral boson masses ]T|| 

The remarkable experimental success[] of this relation which has been proved 
by the precision tests at LEP and SLD, makes the above 50(4) symmetry 
unique and unavoidable for any realistic description of the electroweak sym- 
metry breaking. The same success is far from being achieved in the fermion 
sector, because there are too many models and symmetries which roughly 
give acceptable predictions (i.e. in terms of order of magnitudes). The large 
number of fermions and the lack of direct measurements of all the mass ma- 
trix entries make particularly hard to unambiguously select the underlying 
symmetries responsible for fermion mass hierarchies: measuring the CKM 
matrix elements, one fixes only the product of two unitary matrices, coming 
from the separate diagonalizations of the up and down 3x3 quark mass, and 

1 After the inclusion of some expected and calculable radiative corrections. See || and 
references therein. 
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not both of them. 

To cope with these ambiguities, one usually assumes a fundamental symme- 
try and chooses a breaking pattern. The mass matrices obtained are diag- 
onalized. Eigenvalues and eigenvectors are compared with the experimental 
data, i.e. the CKM parameters and the quark masses. This procedure can 
be used to rule out a model or a symmetry if theoretical predictions and ex- 
perimental data disagree. It is, of course, unable to prove that the assumed 
symmetry is the one which is preferred by the data. Ideally one could 
take a completely model independent approach and prove that a symmetry 
is among the best ones that can predict the correct mixings and masses spec- 
trum. This goal is probably unachievable, but a step in this direction can be 
done if one accepts a certain amount of model dependence, selecting a large 
class of models, and then extracting directly from experimental results the 
preferred models within that class. 

In the following, we explicitly discuss how to approach this goal, introduc- 
ing model independent parameterizations for the up and down quark mass 
matrices. The needed number of free parameters will be in general larger 
than the number of experimental constraints and this will be the main lim- 
itation for a direct fit of parameter values^. Instead of minimizing the \ 2 
with respect to the full set of free parameters (which is too large to be man- 
ageable), we get rid of several unknown parameters, treating them as real 
theoretical errors, through a Monte Carlo procedure. 

To exemplify the proposed method, we will focus on two different pa- 
rameterizations inspired by two classes of models: those based on Abelian 
symmetries and those exploiting a horizontal U (2) symmetry. 

2 The experimental constraints 

Values of the six quark masses and of the four CKM matrix element param- 
eters, used in the present analysis are given in Tables [I] and |2|. 



2 The minimum of the x 2 would be not a single point but obscure and intricate surfaces 
in the n-dimensional space of the n parameters. 
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Quark flavour 


Reference value 
GeV/c 2 


m q (m w ) 
GeV/c 2 


top 
bottom 
charm 
strange 


M pole = 1743±5 A 

m b (m b ) = 4.19 ± 0.05 
m c (m c ) = 1.304 ±0.027 
m s (m T ) = 0.125 ±0.030 


175.3 ±5.1 ±0.2 
2.906 ±0.035 ±0.031 
0.667 ±0.014 ±0.023 
0.072 ±0.017 ±0.002 


up and down quarks 


Reference values 


m q {m w ) MeV/c 2 




2* = 0.46 ±0.09 
Q 2 = 4m i" ( r + ^ )2 = 22.0 ± 0.6 


md(mw) = 3.7 ± 0.9 
m u (mw) = 1-7 ± 0.5 



Table 1: Values of quark masses used in this analysis. The second column 
gives the reference values from which they are evolved using the MS scheme. 
For quarks heavier than u and d the last column gives the values extrapo- 
lated at the scale of the W boson mass, the first uncertainty corresponds to 
the uncertainty on the reference values and the second uncertainty is due to 
the extrapolation. Uncertainties quoted for u and d-quark flavours are only 
indicative as mass ratios are directly determined. 



CKM parameter 


Reference value 


A 


0.2237 ±0.0033 


|V c6 | 


(41.0 ± 1.6) x 10~ 3 


P 


0.225 ±0.038 


V 


0.317 ±0.041 



Table 2: Values of the CKM matrix parameters used in the present analysis 

m 

2.1 Quark mass determinations 

Values of quark masses span a large domain ranging from a few MeV/c 2 , 
for u and d flavours, to more than 100 GeV/c 2 for the top quark. Current 
quark masses, defined in the MS scheme have been used and their values 
have been given at the scale of the W mass. It can be noted that the present 
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analysis depends on quark mass ratios and that these quantities are almost 
scale independent. The running of quark masses has been evaluated following 
the work of |J. The value adopted for the strong coupling constant, at the 
scale of the Z boson mass is a s (mz) = 0.1181 ± 0.0020 which corresponds to 
A( 5 ) = 209.5+|^ MeV in the MS scheme. 

The top-quark mass has been measured from the direct reconstruction of 
its decay products 0]. It corresponds to the pole mass value: 

M p t ole = (174.3 ±5.1) GeV/c 2 . (3) 

The expression relating m t (Mj° ie ) and M pole can be found in [[J. 

Lattice determinations of the b and c quark masses can be found in |J. 
The 6-quark mass is obtained using QCD sum rules for the T masses and 
electronic widths following the work of || . Recent evaluations and complete 
references can be found in Q giving: 

m b (m b ) = (4.17 ± 0.05) GeV/c 2 . (4) 

A new analysis || of the energy dependence of hadronic production near the 
threshold for 6-quark production gives a similar value: 

m b (m b ) = (4.209 ± 0.050) GeV/c 2 . (5) 

In the following the value m b (m b ) = (4.19 ± 0.05) GeV/c 2 has been used. 

The c-quark mass was also determined using QCD sum rules applied to 
the charmonium system, for a recent review see || in which they obtain: 

m c (m c ) = (1.23 ± 0.09) GeV/c 2 . (6) 

The same analysis || which was applied to the 6-quark mass determination 
was done in the charm threshold region and gives: 

m c (m c ) = (1.304 ± 0.027) GeV/c 2 . (7) 

This value is compatible with Equation (||) and is more precise; it has been 
adopted in the following. 



The value of the strange quark mass has been obtained [|Kj] |]TTJ from 
analyses of r decays into final states with an odd number of kaons which 
give, respectively: 



m s (m T ) = (120 ± ll exp ± 8|v us | ± I9 th ) MeV/c 2 (8) 
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and 

m s (m T ) = (130 ± 27 exp ± 8| Vus | ± 9 th ) MeV/c 2 . (9) 

The different values quoted for systematic uncertainties of theoretical origin 
depend on the order of the moments of the spectral function which have 
been used in the two analyses. In the following the value m s (m T ) = (125 ± 
30) MeV/c has been used. It corresponds to a central value of 159 MeV/c 2 
when evaluated at the scale of 1 GeV. 

Finally values for light quark masses can be obtained from chiral per- 



turbation theory [O]. Recent developments have been included [O in the 



values adopted in the following: 

^ = 0.46 ± 0.09; Q 2 = ~ K + m d) 2 = 22 Q ± q g (1Q) 

From these values, the ratio between the masses of the s and d quarks is 
obtained:^ = 19.5 ± 1.2. Using the determination of m s {mw) given in 
Table |T], values are obtained for light quark masses at the scale of the W 
mass. 



2.2 CKM matrix parameters 



In the Wolfenstein parametrisation of the CKM matrix ||14|| , the four param- 
eters are designed, usually, as: A, A, p and r\. In the following, in place of A, 
the quantity, |V C &| = AX 2 has been used instead to avoid correlations. 

Updated determinations of their values, which are given in Table 0, can 
be found in fll~5" |. The values of p and i] have been determined within the 



Standard Model framework using the constraints on these parameters coming 
from B° and B° oscillations, b — > u semileptonic decays and from CP violation 
in K and B mesons 0. 



3 This last constraint, which was not available in |iq | has been included in the present 
analysis. 
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2.3 Procedure adopted to apply constraints 

To apply the constraints related to quark masses and CKM elements we have 
defined the function 

(n th — n ex P\ 2 

where 0^ xp ± a^ xp are the experimental inputs as shown in Tables 1 and 
2. The function x 2 {O th ) will be used afterwards: by means of a weight 
exp(— x 2 /2), we will select models with predictions th , close to the experi- 
mental data. Note that for simplicity we have neglected correlations among 
different experimental data. 



3 Looking for model independent parameter- 
izations 

Low energy fermion masses arise from the Yukawa interaction with the light 
Higgs bosons: 

L Yuk = Y* j H u QiUj + YfHiQidj + h.c. (12) 

In the Standard Model the Higgs bosons H d and H u , are the same particle, 
namely Hd = io"2-ff* but, in general (e.g. in super symmetry), they are two 
distinct fields. The Yukawa couplings Y % \ unconstrained and incalculable in 
the Standard Model, arise from more fundamental high energy Lagrangians. 
The Lagrangian given in (|T2| ) is an effective Lagrangian, where only light 
particles appear as physical fields: at such low energy, heavy particles can 
only appear as virtual internal propagators of Feynman diagrams of the full 
theory. For example a process with n light particles ip that goes into k 
light particles through the virtual exchange of one heavy field F (with 
mass m) through the interaction g\ip n F + g2<fi k F can be described by one 
single operator gig 2 ip n <j) k /M 2 containing only the light particles. Similarly, 
the Yukawa interactions in equation ( fLlp can arise from higher dimensional 
operators H u QiUj(p k /M k where the field acquires a vev <fi and thus = 
(p k /M k . In general the heavy particles content could be very complex, but if 
the original Lagrangian is invariant under some (flavor) symmetry, the above 
mentioned low energy operators must obey the same symmetry. Having in 
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mind this physical mechanism, we will show in the next section how some 
specific symmetries can lead to the observed fermion mass hierarchies. 



3.1 Abelian symmetries 

In the following we will introduce a parameterization for the quark mass 
matrices which is inspired by models with Abelian symmetries. To exemplify 
the physical origin of such parameterization we discuss the simplest case of 
one additional U(l)x at the (very high) scale A new ; after we will comment 
on the more general case. 

This extra U(l)x must be broken, and thus we need at least one Higgs 
boson h with i7(l)x-charge q 7^ 0. This field h must also be a S77(2) wea k X 
U(l) singlet, otherwise the electroweak group would be broken at too high 
energies. In the quark sector we call qf L , qf R , qf L and qf R the [/(l)x-charges 
of, respectively, the up and down type quarks (left and right) of the i-th 
family. Also the light Higgs bosons (responsible for the electroweak symmetry 
breaking) can have non-zero charges q Hu < d with respect^ to the new U(l)x 
symmetry. At the very high energy scale, the Lagrangian includes some 
higher order effective operators^ 

^Yukawa = 9$H u Q iUj h n « + tij 1 1 <lQ ,<> jl ' + h.c. (13) 

The constants g l i would be calculable if masses and interactions of the full 
theory were known, otherwise we can only estimate their order of magnitude 
by simple dimensional analysis 

9u >d - (14) 

where M is the mass scale of heavy particles and the exponent is such 
that the mass dimension of the Lagrangian ( |l"3"l) is equal to four (the field h 
has mass dimension 1). In equation (|T3"| ) we have listed all possible operators 
that are responsible for the generation of quark masses. Since the field h is 

4 As already mentioned, in the Standard Model, we have q H " = —q Hd because there is 
only one Higgs field. 

5 Here we have not included operators containing both h and charged conjugated fields 
h* . This is mandatory in supersymmetric theories, in non-supersymmetric models these 
terms will be considered as subleading. 
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a SU(3) co \ OT x S77(2) wea k X U(l) singlet, these operators are invariant with 
respect to the ordinary gauge transformations for any value of (nfj, nfj), but 
if we also require the U(l)x invariance only few values of nfj, nfj are allowed; 
namely those which satisfy the equations 

q Hu -q- L + q] R + < 3 q = o (15) 

q H *-qt L + qf R + nt jq = (16) 

in order to guarantee that the sum of field charges is zero for each Yukawa 
interaction in (|13D. For simplicity, in the following, we will consider only 
the case in which all solutions of the equations ( |T5| , |i6D correspond to non- 



negative integers^, n%j > 0. Solving equations (|15| , [T6|) , one can build an 
invariant operator Qqh nii . In all other cases no invariant exists. 

The breaking of the U(l)x symmetry is expected to occur at very high 
energy, where the field h acquires a vacuum expectation value (vev), i.e. 
(h) = vx 3> Mw ; while it is at the weak scale that light Higgs(es) take the 
vev's (H u ) = v u , (Ha) = Vd respectively. As a consequence, the Lagrangian 
( |13D and the estimate (|T4]) yield the following masses for the up sector 

~M 

and a similar expression for the down sector. We remind that qn and qjR are 
respectively the left-handed and right-handed fermions ?7(l)x-charges, which 
are in general different — q^ ^ 0). If e = vx/M <C 1, then equation 
( \T7\) implies a definite hierarchical structure of the entries in the mass matrix, 
whose orders of magnitude are proportional to the difference between the 
quark [/(l)x-charges q^ L — q™ R . When all charges are known, applying (|17D , 
one predicts the order of magnitude of all masses and mixings. The converse 
statement is not true: even if all mass matrix entries are known, only some 



m ij = 9u v u (v x ) « ^ v u e 5 ; with e = — (17) 



linear combination of the g" L R can be derived from (]17|) . It is preferable to 
replace the qf L R by the following parameters 

U = e^L-1lR+l3R-1iL)/l^ C = e^L-l2R+l3R-l3L)/<l 

(18) 

+ — 1:iL 1lR U _ 13L 12R 

U n u _ n u _|_ nU —a u U n u —n u A-n u —n u 



6 If riij < one can build the invariant operator Qq(h) n%3 where h is the charge con- 
jugated of h. However, in supersymmetry this operator is not allowed for the holomorphy 
of the superpotential; thus, the corresponding entry in the quark mass matrix would 
be zero. 
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With these definitions, the up quark mass matrix becomes 



iSl~iSr+1h 



a\ u tu c 1 ~ ku 



ale 
< c ku 



o 6 c 



(19) 



A similar matrix can be written for the down quarks, with four additional 
parameters defined from the down quark analogue of equation (|T8|) 



a\ d 



af d td s 
a% d td 



i-fcd 



ii d}- td s kd 



a 5 s 
ai s kd 




(20) 



For each matrix entry in ( |H||2"0"| ) , a complex coefficient a^' a with |a"'"| — 1 
has been included to take into account theoretical uncertainties coming from 
the very heavy particle states (see equation (|T7[), and the discussion above), 
which are intrinsic of any low energy effective description. The matrix in 
equation (^) has three different eigenvalues; it is easy to prove that their 
order of magnitudes are equal to the three diagonal elements in fllTf). The 
variables u and c are consequently constrained to be of the order of the 
ratios m up /mt op and rn c harmlfntop respectively. The off-diagonal elements in 
equation ([T9| ) cannot be estimated by simple arguments; in fact deriving the 
exponents t u and k u in equation (^) requires a more complex and rigorous 
analysis and this will be done in the next section. 

The parameterization ( |T9"D has been inspired by a model with just one 
additional U(l)x symmetry and one heavy Higgs h\ its extent of validity is 
much wider but it does not cover the full set of models based on Abelian 
symmetries [0, 17[ |18, |19, 20| . To simplify our analysis, we will restrict 



ourselves to the class of models defined by equation (|19D , which hereafter 
we will simply call Abelian symmetries. A physically interesting example 

with three additional U(l) 



is, m 



is provided by the model studied in 
factors and three new Higgses. Two of these U(l) can be written as a linear 
combination of the £7(l)'s in the Cartan sub-algebra of Eq D SU (5) x U(l)i x 
f/(l)2, and they distinguish families. The third U(l)s is anomalous [|T^, |TTf |. 
The quark matrices that immediately descend from this symmetry choice are 



A 8 A 5 A 3 
A 7 A 4 A 2 
A 5 A 2 1 



(21) 
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and 

/A 4 A 3 A 3 \ 
M d ~ A 3 A 2 A 2 (22) 
\\ 1 l) 

where A is an unknown parameter. One can check that these matrices agree 



with the parameterizations given in (|19| , |20D , once we fix 



^ = 5/8, k u = l/2, t d = 1/4, k d = 
u = A 8 , c = A 4 , d = \\ s = X 2 . 



(23) 



Other authors |2(J simply add to the usual SUSY SU(5), a U(l) flavour 
symmetry that leads to a better gauge coupling unification and neutrino 
phenomenology. They get 



^ = 1/2, fc u = l/2, t d = 2/5, k d = 0, 
u = A 6 , c = A 4 , rf = A 5 , s = A 2 . 

3.1.1 Fitting the mass hierarchies in the Abelian case 



(24) 



The Method 

The goal is to extract the values of t U}d and k Utd from the experimental mea- 
surements. A direct fit of the data is not possible, since the number of free 
parameters (~ 36) in (|19|,^(]) is much larger than the number of observables, 
six mass eigenvalues plus four CKM parameters. 

The main obstacle comes from the coefficients a"' d , whose phases (and 
size) are not theoretically known. To cope with them, we will treat this 
uncertainty as a theoretical systematic error. Namely, we have assigned a 
flat probability to all the coefficients a^' d with 

| < \al d \ < 2, < arg(a"' d ) < 2tt, 

(25) 

< t < 1, < k < 1. 

u, c, s and d are less than one and we randomly take them with a flat distri- 
bution in logarithmic scale. t U}d and k U}d must satisfy the above constraints 
since (by definition) we choose the entry (3,3) of the matrices in (|l9|,20) to 



be the largest one. For any random choice of the coefficients a"' d , of the ex- 
ponents t U)d , k Uid and of the variables u, c, d, s we get two numerical matrices 
for, respectively, the up and the down sectors. The diagonalization of these 
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matrices gives six eigenvalues, corresponding to the physical quark masses, 
and two numerical unitary matrices whose multiplication yields the CKM 
matrix. We have collected a large statistical sample of events. Each one 
of these events can be compared with the experimental data (see section 2) 
through the x 2 '- the event is accepted with probability : 

P(Cf ) = e ~h 2 (of) (26 ) 

where the \ 2 is defined in (|Tl~|). Before applying the experimental constraints, 
events are homogeneously distributed in the variables t Ut a, k Ut d, and proba- 
bility distributions are flat; but after, applying the weight corresponding to 
eq. (p6|) , only points lying in well defined regions of the space t U) d, k Utd have 
a good chance to survive. 

Let us better clarify the reason for such not uniform distributions and the 
physical interpretation of the density of points per unit area in Figures 1-2,4- 
7. Let us assume two different choices,[] of the exponents t u ^ and k Uj d, that 
we call model 1 and model 2, lying in two different regions (1 and 2) of the 
k u ,d vs t U} d plane (fig.l); the Monte Carlo generates two samples of up/down 
matrices, through equations (|TP|,|2D|). Only a fraction p\ (and p 2 ) of matrices 
of the sample 1 (and 2) will pass the experimental constraints (that is eq. 
(|26|)): pi (p 2 ) is the probability that model 1 (2) predicts masses and mixings 
in a range compatible with the experiments. Then pi and p 2 are, respectively, 
proportional to the density of points in regions 1 and 2. From them we can 
argued that the model 1 is P\jp 2 more (or less, if Pi/p 2 < 1 ) likely than 
model 2. Even if our Monte Carlo approach favours most predictive and ac- 
curate models, we also emphasize that one should not mistake these results 
with true experimental measurements. They only give us "natural" range of 
values for the exponents t u ^ and k u ^- 

If Pmax is the density of points at the maximum, we call R = p/p ma x the ratio 
of the probability with respect to the value at the maximum. By definition 
R < 1 . We will show regions corresponding to different values of R. 

The results 

Figures 1 and 2 show the density of points in the planes kd vs td and k u 
vs t u . Points are blue (black) in regions where R > 0.1 and green (light gray) 

7 We will often use the word "model" to understand a particular and fixed choice of the 
exponents t u> d and k Ut d- 
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1 



down vs t-down 



0.8 




Figure 1: Abelian Symmetries. The exponent k d (vertical axis) vs t d of the 
parametrization given in (POP- Points are blue (black) (green (light gray)) in 
regions where R > 0.1 (R < 0.05) respectively 



where R < 0.05. 

A large fraction of the t and k space is essentially forbidden. We can put 
an upper bound to all the exponents. Namely, t d < 0.55, k d < 0.4 t u < 0.7, 
k d < 0.65. This is expected since large values for t u ^ d and k Utd increase the 
off-diagonal elements in the CKM matrix, in contrast with the experimental 
observation. The model previously mentioned |18|, [19|, prescribing eq. (|23|), 
lives on the border between the highly populated and the empty regions: in 
the up sector (fig. 2) a slightly higher t u would fall in the not populated region. 
Similar conclusions, but in the down sector (fig. 1), apply for the model [p0|] , 
with exponents given in equation (p4|). In these specific models, a closer 
study of the evolution of the Yukawa couplings (through the renormalization 
group equations) from the unification scale down to the weak scale could be 
desirable if one wishes a more definite conclusion. 

One should also note that, even if the experiments put strong bounds on the 
CKM parameters, the allowed regions for t u ^ d and k u ^ d are quite large; we 
have verified that a reduction on the errors of the p and 77 parameters would 
have a negligible impact on Figures 1 and 2. We conclude that errors coming 
from the theoretical uncertainties ( p5|) are dominant. 
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kup vs t U p 



. 8 




0.6 



. 4 



0.2 



0.2 



0.4 



. 6 



0.8 



1 



Figure 2: Abelian Symmetries. The exponent k u (vertical axis) vs t u of the 
parametrization given in flT5p. Points are blue (black) (green (light gray)) in 
region where R > 0.1 (R < 0.05) respectively 

For the moment, a large range of possibilities exists: even a scenario with 
all exponents set to zero t Uj d = k u ^ = like in the following mass matrices 
gives acceptable predictions. 



3.2 U(2) horizontal symmetry 

Another class of models is based on a U(2) horizontal symmetry. This sym- 
metry acts on the known fermion families as follows. The light quarks trans- 
form as doublets 2 under the U (2) group 




(27) 




(28) 
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ft are the lefthanded SU(2) mea k doublets, while f R are the charge conju- 
gated of the righthanded SU {2) wea k singlets. The light Higgses (responsible 
for the electroweak breaking) are singlets as well as the quarks of the third 
generation. The U(2) (differently from SU(2)) includes a £7(1) phase trans- 
formation: the 2 and the 2 are not equivalent representations; in particular 
such a U{2) forbids Yukawa interactions like [fJ|, Q 

g u H u u L u R + g c H u c L c R (29) 

as well as all possible mass or mixing terms concerning the two lightest 
generations. On the contrary the top and bottom quarks can have mass 
since they are singlets under the above U(2). To allow the lighter fermions 
acquiring a mass, we need to break the U(2) symmetry in two steps. Firstly, 
the breaking of U{2) — > £7(1) can be induced by a £7(2)-doublet <p a , the 2, 
and a triplet the 3. Exploiting the £7(2) symmetry we can always rotate 
their vev's in order to obtainP] (0 2 ) = ($22) = v 2 and = ($11) = 
($12) = ($21) = which implies the following Yukawa couplings 

gH u c L t R <p 2 + g H u c R t L (f)* 2 + g" H u c L c R $ 22 + h.c. (30) 

or in terms of the quark mass matrix 

/0 \ 

M up cx 1/M v 2 v, . (31) 

An analogous matrix arises for the down sector. At lower energy also the 
£7(1) can be broken by a t/(2)-singlet A a b, antisymmetric under the exchange 
of the indices a and b. This changes the matrix (|3T| ) into 

/0 -v 3 \ 
M=\v 3 v 2 vi (32) 
V v x M u ) 

where the scale of the U(l) symmetry breaking is much smaller than the 
£7(2) symmetry breaking, i.e. v% v 1 ~ v 2 - The zeros in the entries 

M 13 = M31 = M u = (33) 

8 We label the two lightest families with 1 and 2. 
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are a generic consequence of this class of models^ we will exploit the condi- 
tions (|33|) to parameterize the quark mass matrices (we make no assumption 
on the other entries); starting from the texture (|3"3"D, we will considerably 
simplify the problem of extracting all mass hierarchies from the data; and 
at the same time, this will leave us with a reasonably large and assorted 
selection of models. To be more concrete, after the conditions (|33| ) we are 
left with 6+6 non-zero entries that can be parametrized by 12 free complex 
variables as follows (to simplify the notation, we omit the up/down subscript 
in the exponents) 

a%e\~ Pu 
M up = M u I a«ef a$e 2 I , (34) 



, U 



5 C 2 



and 

M down = M d I a d 2 ef afe 4 a\e\ d I . (35) 
V afef 4 ) 

Clearly one could choose a different reasonable parameterization with differ- 
ent parameters: nevertheless the needed 12 parameters would be in one to one 






"l e 3 





d,Pd 
2 fc 3 
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a d 4 € r 4 d 





n d d d 

a 5 e 4 


4 



correspondence with ours (shown in (0,35)) through well defined equations. 
Thus we do not loose in generality, choosing the above parameterization, 
here the only assumption^ is equation (|33"D. The determinant of (0) gives 
the product of the three eigenvalues and thus m up m charm m top ~ eiM^. We 
also observe that M u ~ m top , and e 2 ~ rricharm/mtop (if t + d > 1), then we 
also get that ei ~ m up m charm /mf op . A similar estimate holds for the down 
sector, where e 3 and e 4 are related to the down quark masses. The expo- 
nents p u ,di r u,d, du,d need a more complex analysis that will be done in the 
next section. 



9 Renormalization group effect from the high energy down to the low energy can modify 
this texture. A scenario with (1,3) entries slightly different from zero has been discussed 

by p. 

1( TTowever this strict equivalence will become only approximately true, after the im- 
plementation of the Monte Carlo procedure in the next section. In fact a different pa- 
rameterization would correspond to a different (and non-flat) a-priori distribution for the 
exponents p Ui d, f Ui d, d u ^ (see the beginning of the next section). This ambiguity is the 
necessary prize to pay for our Bayesian approach. 
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The coefficients a%' d take into account the intrinsic theoretical uncertain- 
ties due to unknown fundamental parameters of the high energy physics. 
They are expected to be of order one : |a"' d | ~ 1 with < arg(a"' d ) < 2n. 

The U(2) symmetry breaking, naively discussed above, implies definite 
values for the exponents 

Pu = 1/2, T u = 1, d u = 1 /Q(5\ 

p d = l/2, r d = l, d d = V, { ' 

however a more sophisticated theoretical study could lead to different sce- 
narios, as for instance if one embeds the above picture into Grand Unified 
Theories. Our aim is to extract these exponents directly from the exper- 
imentally measured quantities, in the same line as shown in section (3.1). 
Before moving to this determination, we would like to show with an example 
the importance of the precise determination of the p and fj parameters in 
excluding the above solution flSoT) for the exponents. 



3.3 Importance of the precise determination of p and 
rj parameters 

In this example we use the values for the exponents that correspond to the 
U(2) symmetry breaking given in equation fl3~6|) . This choice implies that 

Hi : 



V cb 
Vts 




(37) 



These relations are not exact due to the presence of the order-1 coefficients. 
On the other hand, as discussed in [fM|, the larger correction is on the first 



relation in (|37j) and is of order 10%. Equation (^) can be used to define 
an allowed region in the (p-rj) plane as shown in Figure (3). This region 
has been obtained by adding an extra 10% random correction to the first 
equation in (|37D to take into account the above-mentioned effect [p4| . In 



this Figure the ~p-rj allowed region is compared with the one selected by the 
measurements of \ek\, |V u &/V C 6|, Am^ and from the limit on Am s . The 
two selected regions are far to be compatible (to be more quantitative see 
next paragraph). The significant improvement in the determination of the 
p and fj parameters, obtained recently, allows to draw clear conclusions on 
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the U(2) horizontal symmetries. As it was anticipated in the high limit 
on the B s -B s mixing parameter, Am s , disfavours the naive U(2) symmetry 
breaking. 

IP" 1 

0.8 
0.6 
0.4 
0.2 

°-l -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

P 

Figure 3: The allowed region for p and fj (the contours at 68 % and 95 
% probability are shown). Coloured contours : regions obtained using the 
constraints given by the measurements of \V u b \ / \ V c b\, \ek\, Am^ and Am s . 
Empty contours : regions obtained using the relation given in equation (|37|) . 

In the next section we extract the exponents of a generic U(2) horizontal 
symmetry directly from the experimentally measured quantities. 

3.4 Fitting the mass hierarchies in the horizontal U(2) 
models 

The a u ' d ~ 1 coefficients in eq. (|3~3| , |3li)D (differently from the Abelian case) 
are redundant: for example, if \d%\ ^ 1, one can reabsorb the norm \a%\ in 
eq. (P^) into a redefinition of r u — > r u + log | a\ \ / log e 2 of the exponent r u ; 
thus one could set \a^' d \ = 1 without loosing in generality. 
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from quark mass relations 
IVubl/IVcbl = ^rn;/m c 
IVtdl/IVtsI = v^/m; 



from £ k , Am d , Am,/Am d 
and IVubl/IVcbl 




Pdown vs Pup 



0.56 



0.46 



0.48 



0.52 



. 54 



0.5 




0.3 0.35 0.4 0.45 0.5 0.55 

Figure 4: U(2) Horizontal Symmetries. The exponent p d (vertical axis) vs 
p u . The letter " U" marks the naive U(2) prediction discussed in the text. 



As in the Abelian case, the range of variation of t u> d is fixed by the require- 
ment that the entries (1,2) and (2,1) are smaller than the entry (3,3) in 
matrices (^,^). d u>d and r u d could go from zero to infinity, but we have 
checked (see after) that for values greater than five the corresponding entries 
were absolutely negligible. In practice, the case r u ^ > 5 (or d u ^ > 5) is 
equivalent to r u ^ = 5 (or d Uyd = 5). 

For any random choice of the coefficients a"' , of the exponents p Ut d, r Ujd , d Ut d 
andQ the vev e's we get six physical quark masses and the CKM matrix. 
The large statistical sample of collected events has been compared with the 
experimental data (see section 2) and a weight defined for each event, cor- 
responding to the probability distribution ( p6|) . We have found that the 

n Note that e\ and £2 in the up sector are independent from £3 and £4 in the down 
sector. All of them are randomly chosen, with flat probability distribution in logarithmic 
scale. 



Namely, we have assigned a flat probability to all the coefficients a"' d with 




19 



^down vs ^"up 




0.5 1 1.5 2 2.5 3 



Figure 5: U(2) Horizontal Symmetries. The exponent (vertical axis) vs r u 

overall probability distribution is maximal in regions shown in Table 3, even 
if it is still high and acceptable in a larger domain. 



Pu 


~ 0.33 


Pd 


~ 0.5 


d u 


>0.3 


d d 


~ 0.15^0.35 


r u 


~ 0.5 




> 1 



Pu 


~ 0.33 


Pd 


~ 0.5 


d u 


>o 


d d 


~ 0.15^0.35 


r u 


>1- 


Td 


~ 1 



Pu 


~ 0.52 


Pd 


~ 0.48 


du 


>0.5 


d d 


0.6 -r 1.2 


r u 


~ 0.5 


Td 


> 1 



Pu 


~ 0.52 


Pd 


~ 0.48 


d u 


>0.3 


d d 


0.6 4- 0.9 


r u 


>1 


Td 


0.6^0.9 



Table 3: Four regions in the space of the exponents p Uy d, r u ^, d u ^ that 
maximize the probability distribution. 

The two solutions on top of Table 3, prefer asymmetric quark mass ma- 
trices: p u ~ 0.33 implies that the entry (1, 2) (in the up sector) is roughly 
the square of the entry (2, 1); also dd < gives an entry (3,2) much larger 
than the (2,3) (in the down sector). 

The remaining two solutions (bottom row of Table 3) appears more sym- 
metric: a slight asymmetry is due to the exponents p u = 0.52 and pd = 0.48. 
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ddown vs Pdo 



wn 



3 




2.5 
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1 . 5 



1 



0.5 



0.425 0.45 0.475 0.5 0.525 0.55 0.575 

Figure 6: U(2) Horizontal Symmetries. The exponent da (vertical axis) vs 
Pd- The letter " U" marks the naive U(2) prediction discussed in the text. 

But these are very close to 1/2, thus giving approximately symmetric matri- 
ces. 

The six exponents show a significant statistical correlation; unfortunately 
we can only project the sample of events on two dimensional plots, thus 
in Figures 4-7, the number of points per unit area tells us the probability 
distribution as function of each pair of the variables p u ^ r U} d, d U: d (while 
the remaining four variables are integrated out); In the Figures , points are 
blue (black) in regions where R > 0.1, red (gray) in regions where 0.05 < 
R < 0.1, and green (light gray) where R < 0.05. In Figure 4 we see the points 
distribution in the plane pd vs p u . Note that the allowed region is small, in 
particular pd is very well determined and R > 0.05 only if 0.46 < pd < 0.51. 
The exponent in the down sector pd is better measured than the up one, and 
its value is close to 1/2 which yields the well known relation 



The parameterization for the up and the down sector are completely equiva- 
lent but experimental data prefer to constrain the down sector, while p u can 
vary from 0.32 to 0.55. With U we have marked the special case of a naive 
U{2) symmetry as in (|36|). 




m s ~ sin# c . 



(39) 
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ddown vs Pup 




0.25 0.3 0.35 0.4 0.45 0.5 0.55 



Figure 7: U(2) Horizontal Symmetries. The exponent dd (vertical axis) vs p u 



In Figure 5 and r u are strongly correlated: either ~ 1 or r u ~ 0.5. In 
fact if Td ~ 1 and r u ^> 0.5 one would have V c b ~ m s /mf,; instead r n ~ 0.5 
and r^>l yield ~ a/ m c /m t . Both are possible and at least one of them 
must hold. In Figure 6 (dd vs pd) one clearly sees that the point U (referring 
to the model (0)) is outside the populated region. One can overcome this 
problem decreasing dd towards the region with maximal density dd < 0.5. 
Note, however, that r d > 1 (fig-5) and d d < 0.5 imply an asymmetric 
down quark matrix (|35|) . In addition, decreasing dd one also pushes pd very 
close to 0.5, and (from Figure 4) this solution slightly prefers p u < 0.4, that 
is the region with maximal density, as also shown in Table 3. Also Figure 7 
put in evidence a similar behavior: p u can easily range from 0.32 to 0.55, for 
low values of dd', instead points with dd > 0.5 live in a much more restricted 
area with 0.5 < p u < 0.55. Before concluding this section it is worthwhile 
to note that our approach can be used to test if a zero entry is favoured or 
disfavoured by data. For example, in Figure 6 one can see that for rf^ > 1 
the shape of the point distribution no longer depends from dd and remains 
constant even for dd going to infinity. This clearly indicates that when the 
entry M 32 = ef becomes negligible (i.e. zero) few acceptable models still 



survive, thus M^ wn = is not excluded by data (provided that pd is slightly 



reduced). 
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4 Conclusions 



The aim of particle physics is to describe all existing phenomena with the 
smallest number of free parameters. In fact one usually expects that a deeper 
understanding of the foundations of the theory automatically implies a reduc- 
tion of the number of free parameters. For the same reason, but in a different 
context, we await for grand unification of forces, which replaces three gauge 
couplings with just one gauge coupling of the embedding (simple) group. 

Differently from the electroweak boson masses, in the fermion sector there 
are several free parameters. Understanding their origin in terms of new bro- 
ken symmetries is difficult because the number of measurements is much less 
than the number of such parameters. The consequent ambiguity in extracting 
the eighteen complex entries of the quark mass matrices makes very difficult 
to disentangle the true underlying (broken) symmetry of the fermion sector. 
An approximate analytical approach, like those based on estimates of the 
observables in terms of powers of the Cabibbo angle, can become inaccurate 
when several physical observables are involved: if a fine tuning of the coeffi- 
cients a's of order 1/3 is required to get the right V c b , one can safely conclude 
that the model is acceptable; but if this tuning has to be repeated for each 
observable, the model becomes unacceptable. We have concluded that a 
global analysis based on a Monte Carlo procedure with model independent 
parameterizations can give us more clear and trustable conclusions. We have 
shown how such an approach can be successfully implemented, achieving 
more information with less theoretical prejudice. The progress brought by 
this Bayesian approach is the fact that we can put some confidence levels to 
some important parameters which otherwise would be unknown. We can test 
how precise are the determinations of the different exponents, and this makes 
clear to model builders which one of the theoretical ingredients is essential 
and which is not. 

In case of Abelian symmetries we have found that the exponents t and 
k are poorly constrained. Still, an upper bound can be given to the four 
exponents: t d < 0.55, k d < 0.4 t u < 0.7, k d < 0.65. On the contrary, no 
lower bound arises from our analysis. The models explicitly discussed in the 
text lie in the highly populated region, even if near its boundary. 
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In the non Abelian scenario, namely when we (only) assume the (1,1), 
(1,3) and (3,1) entries to be zero, much stronger constraints can be set; also 
strong correlations between exponents of the down and up sectors arise. The 
naive U(2) predictions appear ruled out. Instead, two alternative scenarios 
are suggested by current data: one where is preferably around 0.25. This 
value contradicts the naive U(2) prediction d& = 1. At the same time > 1, 
and that requires an asymmetric texture Mg 2 3> M$ 3 - This agrees with the 
analysis p3| . Such large M| 2 , in GUT theories |p6fl , implies a large (2,3) 
entry in the neutrino Dirac mass which, in turn, favours a large mixing in 



the atmospheric neutrinos pO 



A second scenario allows us to keep the same hierachies as in the naive 
£7(2) , but with a slight (left-right) asymmetry in the (1,2) (2,1) entries 
(both in the up sector and in the down sector) . Namely | M± 2 1 ~ 2 1 M 21 | and 
|M*|~0.7|M*|. 
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